Assessment of coastal river water quality in Bangladesh: Implications for drinking and irrigation purposes

Saltwater intrusion in the coastal areas of Bangladesh is a prevalent phenomenon. However, it is not conducive to activities such as irrigation, navigation, fish spawning and shelter, and industrial usage. The present study analyzed 45 water samples collected from 15 locations in coastal areas during three seasons: monsoon, pre-monsoon, and post-monsoon. The aim was to comprehend the seasonal variation in physicochemical parameters, including water temperature, pH, electrical conductivity (EC), salinity, total dissolved solids (TDS), hardness, and concentrations of Na+, K+, Mg2+, Ca2+, Fe2+, HCO3-, PO43-, SO42-, and Cl-. Additionally, parameters essential for agriculture, such as soluble sodium percentage (SSP), sodium absorption ratio (SAR), magnesium absorption ratio (MAR), residual sodium carbonate (RSC), Kelly’s ratio (KR), and permeability index (PI), were examined. Their respective values were found to be 63%, 16.83 mg/L, 34.92 mg/L, 145.44 mg/L, 1.28 mg/L, and 89.29%. The integrated water quality index was determined using entropy theory and principal component analysis (PCA). The resulting entropy water quality index (EWQI) and SAR of 49.56% and 63%, respectively, indicated that the samples are suitable for drinking but unsuitable for irrigation. These findings can assist policymakers in implementing the Bangladesh Deltaplan-2100, focusing on sustainable land management, fish cultivation, agricultural production, environmental preservation, water resource management, and environmental protection in the deltaic areas of Bangladesh. This research contributes to a deeper understanding of seasonal variations in the hydrochemistry and water quality of coastal rivers, aiding in the comprehension of salinity intrusion origins, mechanisms, and causes.


Introduction
Salinity refers to the saltiness or amount of salt dissolved in a water body.These salts include sodium chloride (NaCl), magnesium sulfate (MgSO 4 ), potassium nitrate (KNO 3 ), and sodium bicarbonate (NaHCO 3 ), which dissolve in water to produce ions.The percentage of salt in various waters is categorized as freshwater (less than 0.5%), brackish (between 0.5 and 30%), saline (between 30 and 50%), and brine (more than 50%) [1].Globally, salinity intrusion is a pressing environmental issue, affecting approximately 20 million people through soil and water in Bangladesh alone.The coastal belt of Bangladesh, covering 32% of the country and encompassing 19 districts [2], has seen a 26% increase in salinity over the last 35 years [3,4].The coast's connection to major rivers through various water inlets and estuaries near the Bay of Bengal makes this bay a primary source of saline water.Cosmic factors, such as the positions and movements of the Earth and sun, indirectly influence environmental elements like temperature and sea-level rise, which can affect salinity.However, these processes occur at a very slow pace.Anthropogenic activities, including deforestation, inadequate water governance, industrial growth, and the extensive use of groundwater for domestic and agricultural purposes, are significant contributors to salinity intrusion.These human-induced actions, coupled with poor infrastructure maintenance and cross-boundary river policies, exacerbate the exchange of cations and anions in the salty aquifer, leading to increased salinity [5].During the full moon, tidal pressures can easily weaken embankments, allowing saline water to enter the polders and causing salinity intrusion.This intrusion is exacerbated by factors such as tidal flow, tidal discharge, wind flow, oceanic pressure, sea-level rise, melting glaciers, and the expansion of oceans [6].Bangladesh's coastal areas are particularly vulnerable to disasters due to their geographic location and the tropical monsoon climate, which brings cyclones from May to November.Cyclones, originating in the Indian Ocean, contribute to increased water salinity in these regions.Between 1584 and 2021, 62 tropical cyclones have struck Bangladesh's coastal belt, resulting in significant damage.Notable cyclones occurred in years such as 1584, 1699, 1767, 1822, 1831, 1847, 1876, 1958, 1960, 1962, 1963, 1965, 1966, 1970, 1974, 1991, 1994, 1995, 1997, 1998, and 2007. Specifically, in 1876, 1970, 1991, and 2007, cyclones reached maximum wind speeds of 222 km/h, 240 km/h, and 260 km/h, with surge heights of 13.6 m, 10.6 m, 8 m, and 12 m, respectively, illustrating the severe impact of these natural events on the region [7].
Global temperatures have increased significantly due to economic growth and greenhouse gas emissions, with an average rise of approximately 0.9˚C to 1.2˚C (equivalent to 1.6˚F to 2.2˚F) above pre-industrial levels [8,9].The Ganges-Brahmaputra-Meghna (GBM) river delta receives substantial freshwater inflow from the melting of Himalayan ice.Monitoring of the Mean Sea Level [10] shows that the global sea level has risen by 210 mm over the past 130 years, with the rate of increase accelerating in recent decades.Predictions suggest that by 2050, sea levels could rise by up to 1 meter along the Bangladesh coast, potentially displacing 5.73 million people [9].In the region of Char Changa, an annual geological uplift of 0.018 mm has been observed [11].Consequently, areas like Hatiya have faced regular tidal flooding and riverbank erosion, with a net sea level rise of 2.94 mm over 20 years.Specific locations such as Hiron Point, Char Changa, and Cox's Bazaar have experienced sea level rises of 4 mm/year, 6 mm/year, and 7.8 mm/year, respectively.Bangladeshi rivers have deposited at least 2.4 billion tonnes of sediment annually for the past 5,000 to 7,000 years.This sediment accumulation has contributed to the formation of new charlands, although it also influences local sea level dynamics on riverbeds [12,13].Tidal ranges along India's western and eastern (central coast) borders vary from 3 meters to 5 meters at the Meghna estuary mouth, where the average elevation above sea level is less than 3 meters [14].In the catchment area, which includes canals, inlets, and rivers, salt levels are higher due to changes in fluvial morphology that decrease freshwater discharge and lead to shoreline waterlogging.The balance between groundwater aquifer recharge and outflow is crucial for preventing saltwater intrusion into freshwater systems.When groundwater aquifers deplete faster than they recharge, saltwater intrusion becomes a risk.The tropical monsoon brings high tides and salty water overflows along Bangladesh's coastline, but moderate to high humidity and temperature result in significant rainfall between May and October, mitigating salinity effects.In 2001, the Huanghe Delta extracted 100 million m 3 of groundwater for irrigation and aquaculture [15].Strategies like those in the Mississippi Delta, which involve irrigating with upstream freshwater to prevent agricultural salinization, and the construction of canals and dikes in the Mekong Delta to supply freshwater irrigation, have paradoxically increased salinity.Brackish water from prawn cultivation has exacerbated soil and aquifer salinization in deltas such as Huanghe [16], Mekong, and GBM [17].Freshwater inflow into deltas, affected by upstream water and land management, leads to salinization from upstream to downstream.Dams and water diversions on major rivers (e.g., GBM, Mekong, Rhine, and Meuse) [18], and consumptive agricultural water use significantly reduce inflow to deltas and distributaries globally.Water diversion within deltas, such as in Venezuela's Orinoco delta where a distributary was dammed for hydroelectricity, poses challenges [19].Additionally, the loss of sediment influx in deltas like the Mississippi has led to wetland subsidence, increased risk of storm surge flooding, and groundwater saline intrusion [20].
Surface water chemistry is essential for providing insights into water types, geochemical processes, and classifications based on hydrochemical benchmarks.By integrating surface water chemistry with geochemical characteristics, we can identify trends, environmental challenges, and enhance our understanding of water sources, geochemical processes, water quality, and susceptibility to contamination.The number and types of salts in water are critical factors affecting irrigation quality, with increased salinity, decreased permeability, and the presence of harmful ions significantly degrading water quality [21,22].Therefore, physicochemical properties are utilized to evaluate irrigation quality [23], employing replicative methods such as those developed by the US Salinity Laboratory [24] and the Wilcox Diagram [25].These methods assess the suitability of water for irrigation purposes.Beyond these, Water Quality Indices (WQIs) are instrumental in evaluating the irrigation quality of water.However, relying solely on one indicator for assessing irrigation water's validity can be misleading and result in poor assessment performance [21,26].Numerous studies have advocated for a water quality index that uses component-weighted scores to provide a more nuanced evaluation [26].WQIs, incorporating factors like Sodium Percentage (Na%), Sodium Adsorption Ratio (SAR), Permeability Index (PI), Kelly's Index (KI), and Residual Sodium Carbonate (RSC), are designed to meet comprehensive monitoring and assessment criteria by consolidating various factors into a singular measure.This aims to classify waters based on their potential uses, chemical composition, and physical properties to guide water allocation effectively.In Limpopo, South Africa [27], demonstrated that the IWQ is an effective tool for monitoring the Luvuvhu Catchment, with parameters such as Electrical Conductivity (EC), Na%, SAR, PI, and RSC providing reliable predictions of irrigation water quality [28,29].Additionally, the adoption of multivariate statistical methods, like Cluster Analysis (CA), which groups physicochemical properties based on the interactions between water's chemical components, has gained popularity.Such geochemical and multivariate analyses are crucial for a comprehensive understanding of water quality and the processes affecting especially contaminated surface waters [30].Principal Component Regression (PCR) is frequently used to predict hyperspectral responses in samples, offering valuable insights into the in situ canopy and broader environmental monitoring efforts [31,32].
Prior research has not extensively evaluated the suitability of coastal waters of the GBM delta for irrigation and drinking purposes.This study aims to identify the spatiotemporal variations in salinity intrusion within the estuarine systems of the GBM delta, offering insights into how these changes impact water quality over time and space.Additionally, this study will determine the entropy-weighted water quality, a method that considers the variability and distribution of water quality parameters, to assess its suitability for drinking and irrigation purposes.By integrating these approaches, the study seeks to provide a comprehensive understanding of water quality dynamics in the GBM delta and inform sustainable water management practices.

Study area 2.1.1 Permission for sample collection.
No permission was required to collect the water samples for this study.The samples were collected from the coastal area adjacent to the Bay of Bengal.As the sampling locations are not under any specific governance or private ownership, there was no requirement to obtain permission from individuals or institutions for sample collection.
2.1.2Geographical information of the sampling area.Covering an area of 47,201 km 2 [33], the coast of Bangladesh (Fig 1) is one of the most complex and dynamic deltas in the world [34].The estuarine system of the GBM delta is classified into three parts: the Western Estuarine System (WES), Eastern Estuarine System (EES), and Central Estuarine System (CES).This classification is based on their hydromorphological characteristics, which are influenced by distinct geological formations and seasonal variations [35].
The three main rivers, Ganges, Brahmaputra, and Meghna, carry approximately 90% of the freshwater discharges from upstream, mixing with the saline water of the sea in the estuaries.The region features relatively flat geography, contrasting with the hilly terrain of the Chittagong region that slopes gently toward the Bay of Bengal [4], resembling an inverse funnel.The active lower Meghna estuary, including rivers like Tertulia, Lohalia, Feni, Halda, Karnafuli, Sangu, and Matamuhuri, drains through the combined flow of the GBM river systems.It is predominantly influenced by the Eastern Estuarine System (EES) region and flows into the northeastern corner of the ocean through this estuary.The Buriswar, Bishkhali, and Baleshwar estuaries, part of the Central Estuarine System (CES) coasts, consist mainly of plain land [35].They are interconnected with the Tetulia estuary and divert GBM river channels.The world's largest mangrove forest, the Sundarbans, located in the CES region, is ecologically very sensitive to salinity intrusion [36].The underlying alluvial aquifer systems, hydraulically connected to the complex network of rivers and channels, are influenced by significant tidal activity and sea-level fluctuations, which impact water availability during low tide (1 m) and high tide (3.1 m) daily.The coastal aquifer's lithology comprises thick deposits of silty clay, swampy peat soil, and coarse to fine-grained sandstone, covered with deltaic sediments and alluvial deposits, with the region's aquifer dynamics strongly affected by active neotectonism.As shown in Fig 1 and Table 1, a total of 45 water samples were collected across the coastal belt areas, located between 92˚19'39.3"Eand 89˚17'34.6"Elongitude and 20˚37'50.6"Nand 21˚40'57.8"Nlatitude.

Hydromorphological and hydrogeological settings
The coastal region of the GMB Rivers forms the Bengal Foredeep Basin, with its eastern part known as the Bengal Flood Plain [37,38].This area is subject to various physical and meteorological phenomena, including storm surges, tropical cyclones, southwest monsoon winds, tides, and tidal flows from the Bay of Bengal.These factors contribute to the mixing of chemical compositions in water through ion exchange, with backwater effects from the sea to the river depicted in Fig 2(A).The Muhuri project, located at the Feni river mouth within the Central Estuarine System (CES), introduces a groundwater-based irrigation scheme that interacts with the saline intrusion affecting coastal belt water and aquifers, impacting agricultural practices such as shrimp farming, which has led to increased soil and water salinity.River systems flow into the ocean in three distinct manners, illustrated in Fig 2(A).High flow occurs when the coastal line (hs) exceeds the depth (hn), leading to erosion and increased flow velocity.In the transition area, the depth remains significant along the channel network, facilitating the observation of an offshore river plume, with mean sea level variations influencing river depth.Conversely, when the shoreline depth (hs) surpasses the typical flow depth (hn), low flow and backwater effects emerge, altering the flow area as indicated by the length of the arrows in the figure.These dynamics cause water surface elevation near the river mouth, especially due to transverse plume spreading [39].Seasonally, from June to October, Chandpur experiences a discharge of 2.5 million cusecs, while the estuary discharges up to 4 million cusecs [40].In contrast, winter water flow reduces to nearly 1/8 of this volume due to the river's slowness and the broad estuary.During the rainy season, water flow can exceed 5 million cusecs, highlighting the significant seasonal variability in the region's hydrology [41].
Salt intrusion occurs when saltwater, which is 2-3% heavier than freshwater, extends over tens of kilometers.In this dynamic, freshwater floats on top of saltwater, leading to a stratification effect.Fig 2 illustrates the process and mechanism behind salt intrusion.The mixing of fresh and saltwater, along with counterpressure from freshwater inflow, serves to reduce salinity further upstream.Salt intrusion influences river discharges, sea levels, and the morphology of estuaries.At the juncture where the river meets the sea, saltwater can penetrate the estuary, as shown in Fig 3(A).Technologies designed to halt salt intrusion operate effectively under normal conditions in river estuaries.However, due to the construction of upstream barrages, the rate of river discharge decreases, enabling saltwater to move inland more rapidly, as depicted in Fig 3(B) and 3(C).In Bangladesh, the discharge rate of the Ganges River at the Hardinge Bridge dropped to 150 m 3 /s in 1995 from a previous 2,000 m 3 /s, as referenced in [33,42].Following the construction of the Farakka Barrage, there was a significant alteration in the hydrology of the Ganges.The Ganges-Brahmaputra-Meghna (GBM) river system could accommodate saltwater navigation if fish ladders and ships are able to pass through sluice gates.

Sample collection, preservation, preparation, and analysis
During the pre-monsoon, monsoon, and post-monsoon seasons of 2020, water samples were collected from various coastal regions of Bangladesh.As indicated in Table 1 and Fig 1, a total of 45 water samples were collected across the coastal belt area.Samples S1 through S15 were collected before the monsoon season, S16 through S30 during the monsoon season, and S31 through S45 after the monsoon season.The sampling locations were selected at the river system's low tide.Before storage, each bottle was accurately labeled with the date, location, and time of collection and then stored in a refrigerator at 4˚C in the dark.The collected samples were analyzed using the Laboratory Quality Management System (LQMS) and standard methods, such as those outlined in the American Public Health Association (APHA) guidelines (2017) [43].Specific physicochemical parameters, including water temperature, total dissolved solids (TDS), turbidity, salinity, and electrical conductivity (EC), were measured.Freshwater samples' parameters were measured on-site with a calibrated portable multimeter (Model: SENSION™156, HACH, USA), while pH and dissolved oxygen (DO) levels were determined using pH meters (Model: HannaHI-255, USA) and DO meters (Model: YK), respectively.Sodium and potassium contents were measured using a flame photometer (Model: PFP7, Jenway, UK).Bicarbonate concentration was determined with 0.2 N HCl, using phenolphthalein as the indicator.Total water hardness was calculated using a complexometric method (TH), expressed in milligrams per liter (mg/L) of calcium carbonate.The hardness was calculated using Eq (1).

Hardness as CaCO
Where, A = Volume of titrant (EDTA) is used in mL and B = Molarity of titrant.An atomic absorption spectrophotometer (Model: AAS240FS, Varian, Australia) was utilized to determine the concentrations of other metallic ions (cations), such as Ca 2+ , Mg 2+ , and Fe 2+ .The concentrations of anions such as Cl -, NO

Quality control for accurate analysis
The study's water samples were analyzed comprehensively in accordance with the ISO/IEC 17025:2017 standards.All calibrated glassware, including pipettes, burettes, volumetric flasks, and beakers, was thoroughly washed with deionized water and then rinsed with 2% (v/v) HNO 3 .An electric digital balance (HR-200, Max 210 g, d = 0.1, A&D Company Limited, Japan; Model N92, D0001), which had been calibrated, was utilized for weighing the samples.Deionized water (Barnstead International, USA, Model-D7071, Resistance 18.2 MO�cm, Conductivity 0.2 μS/cm at room temperature), certified reference material (CRM) reagents (Sigma Aldrich, Germany), and CRM stock standards (concentration 100 ± 4 mg/L) were employed for preparing the samples and standards and for calibrating the instruments.Quality control and assurance measures, including spike recovery, independent standard verification, blank checks, and duplicate analyses, were conducted.The correctness and precision of the analytical data were verified through triplicate analyses using traceable CRM standards, samples, and reagents, and a quality control chart was maintained to monitor these analyses.The concentrations of various metals and anions were determined using atomic absorption spectrophotometry (AAS) and ion chromatography (IC) methods, employing calibration standards that were diluted from stock standards.The detection limits for heavy elements such as Na, Mg, K, Ca, and Fe were 1.0, 0.1, 1.0, 0.1, and 0.2 mg/L, respectively.For metal analysis, the digestion process involved taking a 100 mL water sample, adding 2 mL of concentrated HNO 3 (67-70%), and then evaporating the mixture to reduce the volume to 40 mL.The sample was subsequently filtered using filter paper (Whatman 42 μm) after being diluted to the appropriate volume with distilled water [44][45][46][47].

Agricultural suitability analysis
The quality of water for irrigation purposes is determined by several key factors: Soluble Sodium Percentage (SSP), Sodium Adsorption Ratio (SAR), Kelly's Ratio (KR), Magnesium Adsorption Ratio (MAR), Permeability Index (PI), and Residual Sodium Carbonate (RSC).These parameters are derived using specific formulas, and their results are presented in Table 2.The SSP was calculated using Eq (2), as provided by Todd (1980), with the concentration of all ions expressed in mg/L.The sodium hazard was assessed using the percentage of sodium (Na%) [48].

Assessment and characterization of the water quality
The study includes a comparison of analyzed parameters with the permissible limits for drinking water quality set by various agencies, including the Bureau of Indian Standards (BIS) in 2012 [53], the U.S. Environmental Protection Agency (US EPA) in 2018 [54], the World Health Organization (WHO) in 2011 [55], and The Environment Conservation Rules (ECR) in 1997 [56].Additionally, the entropy-based water quality index (EWQI) was calculated using consecutive footprints to facilitate proper evaluation [57].The eigenvalue matrix (A) is revealed and estimated by Eq (8), where the entropy weight 'i' (i = 1, 2,. .., m) represents the number of river water samples, and 'j' (j = 1, 2,. .., n) represents the number of analyzed parameters. A Eqs (9, 10) are used to evaluate the efficiency and normalization of the analysed parameters.The eigenvalue matrix (A) was resolved into a standard grade matrix (B) for assessing the unique measurement units from various measurement units of the analysed parameters [58].

Basic multivariate statistical analysis
The Piper Trilinear Diagram [59] and the Durov Diagram [60] were utilized to categorize the water class and formations based on the major ions present.To construct these diagrams, Grapher software was employed.Initially, a text file containing the concentrations of anions and cations was prepared.This file was then uploaded into Grapher, where all the data was transferred to an Excel sheet within the software.Subsequently, using this plotted data, both the Piper and Durov diagrams were drawn.Similarly, the Wilcox Diagram, which plots Electrical Conductivity (EC) against Sodium Percentage (Na%), was created using a separate software designed for diagram drawing.This involved taking a text file with the concentrations of cations and anions, separated by commas, uploading it, and then drawing the Wilcox curve to assess the river water quality for irrigation purposes [61,62] , that the coastal river water in the deltaic region of Bangladesh predominantly falls into the NaCl and NaHCO 3 categories, with sodium (Na + ) as the dominant cation and chloride (Cl -) and bicarbonate (HCO 3 -) as the prevalent anions in most river water samples [61].Multivariate statistical analyses were conducted using IBM SPSS version 23 for Windows, including frequency distribution (FD) and Q-Q plots, Pearson's correlation coefficient (PCC), principal component analysis (PCA), and hierarchical cluster analysis (HCA).These analyses were essential to provide a multivariate perspective on the relationships between the calculated parameters.Reliability analysis was employed to quantify the scale of relationships between individual items and their reliability.The internal consistency and group relatedness were estimated using Cronbach's alpha test.Cronbach's alpha reliability coefficients are classified as follows: exceptional (>0.90), good (>0.80),acceptable (>0.70), questionable (>0.60), poor (>0.50), and unacceptable (�0.50).Cronbach's alpha values above 0.81 indicate reliable data and strong correlations among water quality metrics for the analyzed samples.Table 2 lists all Cronbach's alpha values.The suitability of data for PCA or factor analysis was evaluated using the Kaiser-Meyer-Olkin (KMO) test and Bartlett's test of sphericity.The KMO test measures the adequacy of sampling for each variable in the model.A KMO value close to 1.0 suggests high suitability for factor analysis, while a value around 0.50 indicates unsuitability.Values between 0.8 and 0.9 are considered exceptional, 0.7 to 0.8 are good, 0.5 to 0.7 are mediocre, and values above 0.5 are deemed acceptable.In this study, the KMO value was 0.83, indicating that the sampling method is statistically significant for factor analysis.Bartlett's test of sphericity assesses the appropriateness of the data for factor analysis.Table 2 presents the results of Bartlett's test (values less than 0.05).The water samples showed a chi-square (χ2) value of 948.99 with 171 degrees of freedom (D f ).The factor analysis is considered valid since the pvalue (Sig.) is 0.000, which is less than 0.05.

Quality assessment of investigated water
Various organizations, including the Bureau of Indian Standards (BIS) in 2012, the U.S. Environmental Protection Agency (US EPA) in 2018, the World Health Organization (WHO) in 2011, and The Environment Conservation Rules (ECR) in 1997, have established standard permissible limits for drinking water quality parameters.No sampling sites exceeded the standard limit for sodium concentrations, which is set at 200 mg/L.Furthermore, all sample sites, S1 through S45, exceeded the recommended chloride concentration range of 150-250 mg/L.During the pre-monsoon and post-monsoon seasons, the bicarbonate levels at 67% of the sampling sites (30 sites), specifically S1-S15 and S31-S45, surpassed the recommended values.The entropy water quality index (EWQI) was utilized to assess the quality of river water by comparing the significance of each analytical parameter [58].According to the BIS (2012), US EPA (2018), WHO (2011), and ECR (1997) guidelines, the EWQI values for each of the river water quality parameters, such as pH, DO, EC, TDS, total hardness, Na + , K + , Fe 2+ , Ca 2+ , Mg 2+ , Cl -, SO , and HCO 3 -, were calculated.Water samples can be categorized into five groups based on the EWQI values, and they are extremely poor (unacceptable, EWQI � 200), poor (not suitable for drinking, 150 � EWQI < 200), industrial usage, 100 � EWQI < 150), moderate (for domestic, irrigation), suitable for drinking (50 � EWQI <100) and excellent for drinking (EWQI <50).The EWQI value was found in the range between 0.013 and 7.026 with a mean value of 20.87, indicating that the investigated water samples are suitable for drinking (EWQI>50).It also illustrates a homogeneous distribution of water quality for all sampling sites.Evaluated parameters usually influenced the water quality with minimal entropy (ej) and entropy weightage (ωj) values [57].Depending on the EWQI values shown in Table 3, the evaluated water quality parameters followed the order of TDS >Turbidity > DO> Hardness > SO 4 2- Based on the EWQI values, water quality from diverse sampling stations disintegrated as S10 The Na + concentration percentage in the analyzed samples was determined to assess their suitability for agricultural use [62].The observed concentrations of Na + and HCO 3 suggest the presence of alkaline water, which may originate from lithological rock, coastal aquifers, and seawater.When mixed with Cl -ions, this water becomes saline.For water to be suitable for agricultural purposes, the maximum allowable Na + concentration should not exceed 15% [48,62].According to the guidelines set by Wilcox (1948) [62] and Roy et al. (2019) [61], most of the water samples analyzed fell into the doubtful category regarding suitability for irrigation (as depicted in Fig 5), with only two samples classified as unsuitable.

Pearson's correlation coefficient
Pearson's correlation was utilized to assess the relationships between parameters in the samples.Pearson's correlation coefficients (r) are considered significant at the 95% confidence level for the t-test and at the 99% confidence level for the alpha test [63,64].The correlation coefficient (r) was classified as "very strong" (0.80-1.00), "strong" (0.60-0.79), "moderate" (0.40-0.59), "weak" (0.20-0.39), and "extremely weak" (0.00-0.19).multivariate dataset [65].According to the Kaiser criterion, eigenvalues greater than 1 are considered significant and retained, determining the number of principal components (PCs) to be kept [65].This criterion is used to maximize the variance explained by the PCs.Table 4 and Fig 7 detail the loadings of the PCs and their respective variances for the water samples analyzed.Loadings with significance values greater than 0.7 are highlighted to indicate their importance.The aim of rotation in PCA is achieved by revealing many near-zero loadings and a few strong loadings, which clearly distinguish each component.To maximize the variance explained by the principal components (PCs), an effort is made to achieve loadings that are extremely high (positive or negative) or nearly zero, as noted by Davis (1986).For the water samples in this study, all three identified components have eigenvalues greater than 1, as detailed in Table 4 (10.1, 3.1, and 1.7).The scree plot ( Fig 7) demonstrates that these three PCs account for 78.3% of the total variance in the data set under investigation.High eigenvalues signify the importance of each component in representing the data, with the total variances explained by PC1, PC2, and PC3 being 52.95%, 16.3%, and 8.76%, respectively, for water quality data.According to the loading values of 0.94, 0.87, 0.94, 0.86, 0.81, 0.88, and 0.83, with an eigenvalue of 10.1, elements such as hardness, chloride, bicarbonate, sodium, calcium, and magnesium are strongly and positively correlated with PC1.This suggests a significant association of these elements in the water, likely due to the dissolution of calcite or dolomite.PC2, with an eigenvalue of 3.1 and a loading value of 0.87, is strongly associated with iron (Fe) alone, indicating its distinctive influence on the water quality.PC3 does not exhibit significant correlations with any variables, highlighting the unique contributions of PC1 and PC2 to the water quality analysis.These findings underscore the significance of these components in classifying the water quality and suggest saltwater intrusion into the coastal aquifer at the studied site, as indicated by these quality parameters.Principal Component Analysis (PCA) and Hierarchical Cluster Analysis (HCA) were employed to identify and categorize spatial similarities among water samples according to their analogous behavior with respect to evaluated water quality metrics [66].The dendrogram resulting from HCA, which utilizes Ward's method and a measure of similarity for evaluation, can be segmented into three main clusters (see Fig 8A).The contaminants of the seawater source indicated in PC1 and Table 4 as Cluster 1 included the water temperature, total hardness, bicarbonate, turbidity, CO 2 , Mg 2+ , Ca 2+ , PO  4) and comprised DO and K + , suggesting that dissolved oxygen and potassium minerals were mixed in the river water.Cluster 3 of the pH, TDS, E.C., and Fe 2+ indicators shows that iron ore is leaching into the water body.A statistically significant correlation was found among the evaluated parameters in the analyzed water samples.This correlation mirrored the similarities between the variance explained by Principal Component Analysis (PCA) and the variance within the clustering of sampling locations as determined by Cluster Analysis (CA), as illustrated in Fig 8B .The concordance between the PCA and CA results supports the connection among inherent water quality parameters.The observed common and significant positive correlation among these parameters, given their sources of origin, is likely due to geogenic activities, seawater intrusion, and the dissolution of carbonate minerals.Hierarchical Cluster Analysis (HCA) grouped the sampling sites into three clusters based on their similarities: HCA1 comprised sites S1, S31, S2, S32, S4, S34, S23, S24, S25, S22, S13, S15, S45, S20, S21, S43, S14, S44, S3, S33, S27, S37, and S38; HCA2 included sites S16, S18, S19, S30, S29, S28, S17, and S26; HCA3 consisted of sites S5, S35, S39, S36, S6, S41, S9, S11, S7, S40, S8, S10, S12, and S42.
The sodium adsorption ratio (SAR) significantly influences the water infiltration rate, making a low SAR preferable for maintaining soil health.The percentage of sodium (%Na) is critical in evaluating the suitability of groundwater for irrigation purposes because high concentrations of sodium ions in irrigation water can interact with soil particles, reducing their permeability [68].When water contains a high Na + concentration, clay minerals in the soil absorb these ions and release Ca 2+ and Mg 2+ ions.This exchange of Na + in the water for Ca 2+ and Mg 2+ in the soil reduces permeability and, consequently, soil infiltration rates.As a result, soils often become compact and hard during dry conditions, while the movement of air and water through the soil is restricted under wet conditions [69].The SAR index is a measure of the soil's ability to exchange Ca 2+ and Mg 2+ ions for Na + ions in groundwater at ionexchangeable sites.This exchange leads to soil particle dispersion and a decrease in infiltration capacity [70].Despite potential issues with high salinity, irrigation water can still be beneficial for the land.The Kelly's Index (KI) was calculated to determine whether groundwater could be utilized for irrigation purposes [71].
SAR values ranged between 16.87 and 21.18 mg/L, with a mean value of 16.83 mg/L.The SAR values for the pre-monsoon, monsoon, and post-monsoon seasons were 16.86, 18.63, and 15.80, respectively, as listed in Table 5.The SAR values indicate that the water quality is not alarming across all samples.However, according to the Wilcox diagram [49], the water samples from the study area are classified as unsuitable for irrigation purposes (Fig 9).Wilcox [49] also provided a classification for irrigation water based on the percentage of sodium (%Na), as depicted in Fig 9 .The mean soluble sodium percentage (%Na) values for pre-monsoon, monsoon, and post-monsoon were 63.0 mg/L, 61.0 mg/L, 67.0 mg/L, and 60.0 mg/L, respectively.Additionally, the mean residual sodium carbonate (RSC) values for the pre-monsoon, monsoon, and post-monsoon seasons were found to be 89.22 mg/L, 103.46 mg/L, 57.01 mg/L, and 107.31 mg/L, respectively.These results suggest that the water is generally unsuitable for irrigation.
Excessive carbonates and bicarbonates, relative to calcium (Ca 2+ ) and magnesium (Mg 2+ ) ions, can compromise the quality of irrigation water by causing the precipitation of these alkali earth metals, particularly Ca 2+ and Mg 2+ .This precipitation of Ca 2+ and Mg 2+ as carbonate minerals may lead to an increase in both the sodium adsorption ratio (SAR) and the concentration of sodium ions [72].Furthermore, the high residual sodium carbonate (RSC) can result in the dissociation of organic matter, which ultimately leaves a black residue on the soil's surface upon drying.This phenomenon has the potential to deteriorate the soil's physical properties [73,74].The mean Magnesium Adsorption Ratio (MAR) values for pre-monsoon, monsoon, and post-monsoon were 34.92 mg/L, 37.53 mg/L, 30.94 mg/L, and 34.43 mg/L, respectively.According to Table 5, the water is considered safe for irrigation as the MAR values are below 60 mg/L.Conversely, the values for Kelly's Ratio (KR) were found to be 1.28, 1.19, 1.66, and 1.16, respectively, indicating that the water is unsuitable for irrigation when the KR is greater than 1.Furthermore, the Permeability Index (PI) values were 144.45 mg/L, 147.27 mg/L, 133.88 mg/L and 153.7 mg/L, respectively.These PI values categorize the water as Class III, indicating it is unsuitable for irrigation purposes.

Conclusion and recommendations
River water is an essential source of freshwater, crucial for supporting aquatic life, irrigation, industrial uses, environmental conservation, and drinking.Agricultural projects and natural disasters, such as hurricanes, can affect the freshwater availability in downstream areas.Consequently, salinity intrusion has been progressively affecting river water in Bangladesh, especially across the coastal belt.Typically, the salinity level of freshwater is less than 0.5 parts per thousand (ppt), but this study found that salinity levels in the coastal belt of Bangladesh ranged from 1.7 to 2.3 ppt.River salinity was observed to be lower during the monsoon season and higher during the dry season, which spans from March to April.Overall, the water quality in the coastal belt of Bangladesh has not exceeded the permissible values for general use.However, the Entropy Water Quality Index (EWQI) and Sodium Adsorption Ratio (SAR) data suggest that the samples analyzed are suitable for drinking but unsuitable for irrigation.Additionally, the Permeability Index (PI), Kelly's Ratio (KR), and Residual Sodium Carbonate (RSC) values further indicate that the water is unsuitable for irrigation or farming purposes.
To mitigate salinity intrusion in the coastal surface waters, sustainable polder development initiatives should be launched, including the construction of twister technology or cement concrete block (CC block) dams, routine dredging of river sediments, and halting the construction of upstream dams.To address these activities by neighboring countries, a strategic foreign policy and diplomatic engagement are essential to foster a win-win situation.Dialogues should be initiated in international forums such as the UN, ASEAN, NAM, and SAARC to secure water rights and justice and to combat salinity intrusion in the Teesta, Padma, and Jamuna rivers.Rainwater harvesting can serve as an alternative water source or a viable technique to alleviate salinity intrusion across the coastal belt during the rainy season.The use of groundwater for irrigation and shrimp cultivation should be discontinued.Balancing the upstream discharge rate through projects like the China Teesta Barrage is crucial.Introduction of salinity-tolerant crops through genetic modification or other sustainable methods is recommended.Nature-based solutions, such as planting Sonneretia apetala and Nipa fruticans trees, should be adopted because these trees not only reduce salinity but also mitigate climate change, protect against natural hazards like cyclones, and lower temperatures.Additionally, their fruits can provide alternative livelihood options for the local community through the production of jam, jelly, pickles, sauce, tea, and more.
Fig 6 demonstrates both negative and positive correlations among quality parameters and computed variables.A positive correlation indicates that parameters sharing common sources are not significantly influenced by additional sources, reflecting similar environmental behaviours.Conversely, a negative correlation suggests that the spatial distribution of parameters is affected by common causes, yet their environmental behaviors differ.For instance, electrical conductivity (EC) correlates with total dissolved solids (TDS), salinity, pH, bicarbonate, Na, chloride, sulfate, and Mg; chloride correlates with salinity, pH, dissolved oxygen (DO), turbidity, bicarbonate, EC, phosphate, hardness, nitrate, Ca, and Mg; bicarbonate correlates with water temperature, salinity, pH, hardness, phosphate, EC, chloride, Ca, and Mg; and Na correlates with EC and TDS.Strong positive correlations indicate a common origin and similar behavior in the examined environment, while moderate positive correlations suggest a slightly different origin.The relationships between related water quality metrics ranged from very weak to weak, suggesting a modest contribution to overall water quality.

Table 2 . Suitability test of data for factor analysis of water samples. Suitability test of data for factor analysis
https://doi.org/10.1371/journal.pone.0300878.t002